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The frustrated quantum spin-one-half Heisenberg model on the triangular and 
Kagome lattices is mapped onto a single species of fermion carrying statistical flux 
9 = it. The corresponding Chern-Simons gauge theory is analyzed at the Gaussian 
level and found to be massive. This provides a new motivation for the spin-liquid 
Kalmeyer-Laughlin wave function. Good overlap of this wave function with the 
numerical ground state is found for small clusters. 
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The 2D spin-i antiferromagnetic Heisenberg model has attracted a lot of interest over 
the last a few years. It is widely believed that the ground state has long-range Neel order 
on the square lattice. 1 On the triangular and Kagome lattices, the situation is much less 
clear due to the geometric frustration. In the case of triangular lattice, Huse and Elser have 
constructed variational wave functions with long-range order and low energies. 2 Numerical 
evidence was found supporting this scenario. 3 Recently however, Singh and Huse 4 have 
calculated the sublattice magnetization using a series expansion method, which is believed 
to be accurate, and found that the ground state is nearly disordered for the triangular lattice 
and strongly disordered for the Kagome lattice due to the large degeneracy of classical ground 
state configurations. 5 

Recently a two-dimensional extension of the Jordan- Wigner transformation, which es- 
sentially treats hard core-bosons (see below) as fermions with flux tubes, or equivalently, 
fermions coupled to a Chern-Simons gauge field, was developed. 6 The advantage of this ap- 
proach is that the unwanted hard-core condition in the boson picture is taken care of by the 
Pauli principle, but as a price one has to introduce gauge interactions to fix the statistics. 
We applied this method to both triangular and Kagome lattices. At mean field level, the 
flux carried by the fermions is smeared out to form a uniform background magnetic field. If 
we neglect the Ising part of the Hamiltonian (which becomes a nearest-neighbor repulsive 
interaction between fermions after the transformation) for the moment, we have noninter- 
acting fermions moving in a constant magnetic field. Numerical diagonalization yields two 
and six Landau subbands, with large excitation gaps 23 at the Fermi level of 2^/3 J and 1.46 J 
for the triangular and Kagome lattices respectively (J is the coupling between neighboring 
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spins). The Hall conductance is quantized in such cases, in the form a xy = where e is 
the charge carried by the fermions. In the continuum v is just the Landau level filling factor 
v — | , where 9 is the statistics angle which is n in this case. However on a lattice, v is one 
of the TKKN 7 integers, not necessarily equal to |. If one examines Gaussian fluctuations of 
the gauge field about its saddle point (mean field), one can integrate out the fermion degrees 
of freedom (which are quadratic in the action) and expand the effective action for the gauge 
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field about its saddle point to second order to obtain 8 



S[A,] = So + J d 2 xdt[^E 2 (x,t) - ^B 2 (x,t)} 
+ yd 2 xdt(v-^)e, uX A»F l ' x + --; 



where Sq is the mean field action, e and x are the mean field values of the long-wavelength, 
low-frequency dielectric constant and diamagnetic susceptibility respectively. As noted by 
Fradkin, 8 the fluctuation is massless if and only if the Chern-Simons term in the action is 
cancelled, i.e., v — | . This is easy to understand in terms of self-consistent linear response. 
Assume there is a long-wave length, low-frequency fluctuation of the density of the fermions. 
Since the fermions carry flux, there should be a fluctuation of magnetic field in the same 
mode. According to Maxwell's equations, there will be a nonzero line integral of electric 
field around any region T: 



where $ and N are the flux and number of particles in that region respectively. Now look 
at the response of iV to this electric field: 



We see the above equations are consistent if and only if v = |. This is guaranteed in the 
continuum, where Fetter et al did find a gapless collective mode. 19 We have computed the 
TKKN integer for the triangular and Kagome lattices using the method of MacDonald 20 
and found in both cases v — — 1^| = 1. Hence we expect the Chern-Simons field to be 
massive and the quantum XY model is therefore likely to have a gap assuming no broken 
translation symmetry. This gap may be stable against the Ising perturbation causing the 
Heisenberg model to have a spin liquid ground state on these lattices. 

By making analogy to the fractional quantum Hall effect (FQHE), Kalmeyer and Laugh- 
lin (KL) suggested a very interesting spin liquid wave function for the Heisenberg model on 
the triangular lattice. 9,10 It was found to have reasonably good energy (about ten percent 
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higher than the best numerical estimate). We believe the massiveness of the Chern-Simons 
theory demonstrated here provides a more fundamental motivation for quantum Hall physics 
in frustrated spin systems. Further more, within the single-mode approximation 12 an ex- 
citation gap in a 2D system requires a Jastrow-like wave function whose square is a 2D 
one-component plasma in order for the structure factor to vanish at small q: S(q) ~ q 2 . 
SU(2) symmetry uniquely restricts the coefficient of the plasma charge to be that given 
by the KL wave function (m = 2). The spin-| excitations argued by Laughlin 11 similarly 
require m = 2 in the Bose representation. These arguments strongly suggest that a spin 
liquid wave function should be of the Kalmeyer-Laughlin type. 

In the rest part of this paper we first briefly review the KL wave function, and prove 
that it is equivalent to a projected underlying fermion wave function. We apply this new 
wave function to the Kagome net, and calculate its energy using the Monte Carlo method. 
Then we calculate the overlap between this wave function and the exact wave functions on 
small clusters. Finally we summarize and discuss our results. 

The Hamiltonian for antiferromagnetic Heisenberg model is: 

H = J £ Si ■ Sj, (1) 

<ij> 

where J > 0, Si is the spin operator at site i, and the sum is over nearest neighbors. 
Following Ref[9], we map the spin operators to hard-core boson operators. The Hamiltonian 
in this representation is 

H = T + V 

T = J - ( a l a j + a ) a i) 

<ij> 

V = J ^2 riifij + const 

<ij> 

Here N s is the number of sites of the lattice, aj is the hard-core boson operator at jth site. 
As was shown by KL, 9 on the triangular lattice, after an appropriate gauge transformation, 
this new Hamitonian describes hard core bosons moving in a uniform magnetic field with 
field strength one flux quantum per unit cell, under the symmetric gauge. The bosons have a 
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nearest-neighbor repulsive interaction besides the hard core condition. In the ground state, 
the number of bosons N = N s /2, which means the Landau level filling factor is one half. By 
making analogy to the FHQE, they suggested the trial wave function for the bosons: 

*(*!,•••,**) = II(*-*;) 2 II G{z k )e-W%. (2) 

i<j k<N 

Here Zj = Xj + ii/j is the complex coordinate of the lattice site occupied by the jth boson, 
l = b(V3/4n)^ is the magnetic length which is set to 1 afterwards, b is the lattice constant. 
G(zk) = ±1 are gauge phases introduced by the gauge transformation. 9 

The wave function (2) has some nice features, including being a singlet state in the 
thermodynamic limit, but it can not be generalized to non-Bravais lattices, such as the 
Kagome lattice. Also it becomes a singlet only when the system is infinitely large, so it is 
not suitable for finite size studies. For these reasons, we want to find a more general wave 
function that returns to (2) on the triangular lattice in the thermodynamic limit, and has 
better finite-size properties. To do that, we assume the spins are carried by spin-| fermions 
that are trapped on lattice sites, and try to describe state (2) in terms of these underlying 
fermions. 21 We can express the spin operators in terms of these fermion operators: 

a ) = S t = c h c Jl 

a 3 = = C )i C j] 

where ct^ are creation operators of up(down) spin fermions at jth site. Just as for the 
bosons, there is a hard core condition on the fermions: + nj± = 1. 
In the second quantized representation, the state (2) is 

l*>= E *(2i,---,^)4 1 ---4jo b ). 

{zi,---,z N } 

Here the sum is over all possible boson configurations. |0&) is the boson vacuum state. Since 
it corresponds to the state that all spins are down, we have 

|o*> = ]»/> 
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where |0/) is the fermion vacuum state. So we get 

l*> = ^j E •••,**) 

{zi,—,zn;z[i],—,z[n]} 

x 4 1T c, u •••4 JVT c 2iVi n 4j°/>- 

fc=l 

Here = zjy + j denotes the coordinate of jth down spin fermion. The sum is over all 
possible fermion configurations satisfying the hard core condition. We do not distinguish Cj 
and c z , ctj and a z , etc., if z is the complex coordinate of jth site. Rearranging the order of 
fermion operators and neglecting constant factors, we have 

l*>= E *(*!,•••,**) 

{zi,—,zjv;«[i],— ,«[jv]} 

N N 

x F(z ir - -,z N ;z {1] ,- ■ -,z [N] ) J] 4 fcT II 4 u l°/>- 

k=l 1=1 

Here F(zi, • • •, Z[n]) is a totally antisymmetric factor: 

\F(z±, ■ ■ •, z 2 n)\ = const 
• •> • • •, ^fc, • • z 2 n) = —F(- ■ •, z k , ■ ■ •, z,-, • • •). 

We can take F to be 

F = Y[( Zi - zj)- 1 ^- z {j] )- 1 H fa-z^- 1 . 

i<j k,l<N 

If we go back to first quantization, the wave function that describes the underlying fermions 
is just 

*/(*i> • • •> w = n (*< - z j)( z w - z m)~ l 

i<j<N 
N 

x [] (**-*w) -1 II G(z m )e-il^l 2 . (3) 

fc,KJV m=l 

Obviously a fermion wave function should be antisymmetric under exchange (including 
spin variables), but here we neglect spin variables in the wave function and treat up and 
down spin particles as if they were distinguishable, so the above wave function is legal. 18 
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Since we have the hard core condition on fermions, what we really mean by (3) is the state 
that is projected to the subspace with no double occupancy. Hence it is a well defined wave 
function for the original quantum spins. 

Now we can use a theorem proved by KL 10 : 

where £j is the complex coordinate of jth site, Co is a constant, and G(£fc) are the gauge 
phases. 9 This theorem holds only on the triangular lattice in the thermodynamic limit. Using 
it in (3) we get up to a constant factor, 

• • •, z [N] ) = [] ( z i - z j)( z ii] - z m)- ( 4 ) 

i<j<N 

The right hand side of equation(4) is just the product of two Vandermonde 
determinants, 21 which is what one gets when both spin states of the first Landau level 
are fully occupied. Since up and down spin particles occupy the same spatial Slater deter- 
minant, the resulting state must be a spin singlet, even after projection. 16 This provides 
another way to prove that the state (2) is a singlet in the thermodynamic limit. The advan- 
tage of (4) is it gives a singlet even on finite size systems, and it can be generalized directly 
to non-Bravais lattices. 

As a test of the equivalence between (2) and (4), we calculated the energy of (4) on a 
triangular lattice. The calculation of the Ising part of the energy is straightforward using 
Monte Carlo, and the total energy is exactly three times that (for any system size). The 
extrapolated result is —0.48 ± 0.01J per site, which agrees with the KL result. 9 We have 
found that in our case the data converges much faster, i.e., the finite size results are much 
closer to the extrapolated result. This tells us that (4) is better for finite size study. The 
energy we get for the Kagome is —0.399 ± 0.001J per site, about eight percent higher than 
the best numerical estimate. 15 

We have also studied the overlap between (4) and exact ground state on small clusters, 
where we need to minimize the finite size effect by applying periodic boundary conditions 
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(PB). The wave function (4) satisfies open boundary condition, so we need to solve the wave 
functions with PB: 

t(L j )rl>(z) = e i +trl>(zy,j = l,2,..., (5) 

Here t(Lj) is a magnetic translation operator. 13 This problem was solved for the torus 
geometry by Haldane and Rezayi. 13 It turns out that the filled Landau level wave function, 
after neglecting constant factors, can be expressed in terms of elliptic theta functions 13 ' 14 : 

*/= II Pi(f (*-*i)|r)0i(f (z b] - z b] )\r)} 




Here 6i(z\t) is the elliptic theta function, r = L 2 e t(J /L 1 , L\ and L 2 are the vectors that 
determine the shape of the parallelogram, and S is the angle between them. Z is the center 
of mass coordinate which is determined by 4>i an d 02- In most cases we are interested in, 
it should be set to zero. 17 Like eq. (4), eq. (6) also describes a singlet state. A truely 
nondegenerate ground state wave function must be real. 9 So instead of using the complex 
wave function (6) directly, we use ^e'* + \^e -1 ^ in the overlap calculation and use as 
a variational parameter. We applied it to several clusters of the triangular lattice with the 
shape of a parallelogram (torus geometry) and one with the shape of a hexagon (twelve 
spins). The results are listed in table 1. We have found that the square of the overlap 
remains large in systems with up to twenty spins (4x5). The energies one gets are close to the 
Monte Carlo result, which means the change of boundary condition does not change the short 
distance correlations. For reasons we do not understand yet, the overlap is exactly zero in the 
4x4 cluster. We have verified that \& / has the correct symmetries (spin rotation, translation, 
180° rotation and mirror reflection of space). Apparently there is some additional hidden 
symmetry which does not match that of the numerical ground state. We have done the 
same calculations on Kagome clusters 15 with twelve and eighteen spins. Again we got zero 
overlaps, probably for similar reasons. The energies, — 0.420J and — 0.418J per site for 
twelve- and eighteen-spin clusters respectively, are close to the Monte Carlo result. 
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If the ground state of the triangular lattice has three sublattice order, such order is 
suppressed on clusters 3x4, 4x4 and 4x5 due to incommensurability, but is not on 
6x3 and the hexagon with twelve spins. Our data suggests this commensurability is a 
weak effect. Both state (6) and the true ground states have a lot of symmetries (rotation, 
translation, etc.). It could happen that there are so few states of the right symmetries 
available that an arbitrary combination of them will have a decent overlap with the ground 
state. By assuming singlet states are uniformly distributed in the momentum space, we 
find the number of states with the right symmetry is of order 1000 in the case of twenty 
spins, and yet the square of the overlap is rather large: 0.493. The twelve-spin hexagon has 
additional symmetries, so the significance of the remarkably large squared overlap of 0.966 
is unclear. 

Quantum Hall wave functions exhibit hidden ODLRO due to the binding of vortices to 
charges. 24 This corresponds to chiral order in the present problem. 25 The large overlaps we 
obtain suggest chiral order is present in the exact ground states of the clusters. 

The central result of this paper is the demonstration that treating spins as fermions 
carrying flux tubes leads to a massive Chern-Simons theory on frustrated lattices. This pro- 
vides a new and fundamental motivation for quantum Hall types of spin-liquid physics. We 
have developed the formalism needed to compute the overlap between these wave functions 
and the exact ground state and we have obtained significant overlaps for small clusters. It 
would be highly desirable to see these calculations extended to large lattices, although this 
will require considerable numerical effort. 

We acknowledge helpful discussions with Rajiv Singh, V. Elser, A. J. Berlinsky, C. Kallin 
and N. Read. This work was supported by NSF-9113911. 
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TABLES 



TABLE I. overlaps and variational energies on small clusters of the triangular lattice 

cluster | overlap | 2 E v exact energy 

hexagon(12) 0.966 -0.591 -0.6103 

3 x 4 0.821 -0.519 -0.5776 

4 x 4 0.000 -0.459 -0.5347 
6 x 3 0.554 -0.491 -0.5811 
4 x 5 0.493 -0.481 -0.5581 

Energies are in unit J per site; E v is the variational energy. 
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